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ABSTRACT 

Observations of He-like and H-like Ar (Ar XVII and Ar XVIII) lines at 3.949 A 
and 3.733 A respectively with the RESIK X-ray spectrometer on the CORO- 
NAS F spacecraft, together with temperatures and emission measures from the 
two channels of GOES, have been analyzed to obtain the abundance of Ar in 
flare plasmas in the solar corona. The line fluxes per unit emission measure show 
a temperature dependence like that predicted from theory, and lead to spectro- 
scopically determined values for the absolute Ar abundance, A(Ar) = 6.44 ±0.07 
(Ar XVII) and 6.49±0.16 (Ar XVIII) which are in agreement to within uncertain- 
ties. The weighted mean is 6.45 ±0.06, which is between two recent compilations 
of the solar Ar abundance and suggest that the photospheric and coronal abun- 
dances of Ar are very similar. 
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INTRODUCTION 



The re has been much recent interes t in the abundance o f argon in the Sun and nearby 
objects (ILanz et al.l 120071 : iLodderd 120081 ; lAsplund et al.ll2009l ). The abundance of argon in 
the solar photosphere cannot be directly derived as the temperature is much lower than 
that needed to excite neutral or singly-ionized Ar lines, so proxies such as nearby B stars, 
H II regions, planetary nebulae, or solar system sources such as the atmosphere of Jupiter 
must be used. It is presumed that the Sun collapsed out of material that is represented by 
such objects, so the element abundances derived are "protosolar" : allowance must, however, 
be made for a certain amount of heavy-element settling over the Sun's lifetime. There are 
several emission lines of Ar ions in the extreme-ultraviolet and X-ray spectra of the quiet 
or activ e region solar corona or solar flares, and suc h lines have been used to give Ar abun- 
dances (JFeldman fc Widinglll990l ; lYoung et al.lll997l ). Commonly line fl ux ratios are used to 
obtain the ratios of abundances from spectroscopic measurements; thus iFeldman &: Widing 
(119901 ) obtained argon-to-magnesium abundance r atio of 0.15±0.05 fr om measurements with 
the Skylab ultraviolet spectrometer S082A, while lYoung et al.l ( 119971 ) obtained an argon-to- 
calcium abundance ratio in the range 0.55 — 1.1, from both Skylab and optic al measur e ments 
during an eclipse. Using published element abundance sets for Mg and Ca, ILodderd (120081 ) 
in her review has given a grand average of values from various solar determinations as 
A(Ar) = 6.50 ±0.10. (Abundances here are expressed on a logarithmic scale with H = 12 by 
a notation extensively used in recent literature, A(Ay).) The review by lAsplund et al.l (120091 ) 
gives A(A r) = 6.40 ± 0-1^ usi ng a revised O abundance, which updates a previous deter- 
mination (lAsplund et al.ll2005l ) of 6.10 based on a single measurement from solar energetic 
particles. The possibility of fractionation through the "FIP" (first ionization potential) effect 
is recognized, with the abundances of elements having low FIP (< 10 eV) being enhanced 
in the solar wind and corona by factors of ~ 4 and abundances of elements with high FIP 



-> 10 eV) approximately the same as the proxies for the solar photosphere (JFeldman et al. 



19921 : iFludra fc Schmelall999l ): Ar is of interest because its FIP (15.8 eV) is larger than that 
of any other common element in the Sun apart from He and Ne. 

The high-resolutio n X-ray crystal spectro meter RESIK (REntgenovsky Spektrometr s 
Izognutymi Kristalami; ISylwester et al.l (120051 )) on board the CORONAS-F spacecraft has 
observed X-ray lines in the 3.3-6.1 A range emitted by solar coronal plasmas. Line fluxes with 
diagnostic potential including the possibility of deriving element abundances have already 
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been analyzed (Phillips et al 



20 id ISvlwester et al.l 12010 : 



2003; Sv lwester et al.ll2008al Jbl; ISvlwester. Sylwester. &: Phillips 



Phillips et al.l 120101 ). The instrument operated between 2002 



and 2003, and several thousand spectra were collected. It had significant improvements over 
previous spectrometers in that a background formed by fluorescence of the crystal material 
was eliminated for the two short-wavelength channels by judicious settings of the detector 
high-voltages and pulse height discriminators, and minimized to very small levels for the 
remaining two channels. The period from 2002 August to 2003 February was particularly 
appropriate for analysis as the settings were then fully optimized, and spectra were selected 
from this period for the analysis described here. 

Here we use flux values of the Ar XVII resonance (w) line at 3.949 A (actually the flux 
in the 3.94-4.01 A interval minus the flux in a neighboring portion of continuum) to derive 
the absolute abundance of Ar (Ar/H) by comparing with theoretical values taken from the 
CHIANTI database and code (jDere et al.lll997l , 120091 ). As the sensitivity of RESIK is very 
high, the flux measurements have good statistical quality and are a consi derable improvemen t 
over those from the Solar Maximum Mission Flat Crystal Spectrometer ([Phillips et al.lll994l ). 
The Ly-a lines of Ar XVIII form a single prominent feature at 3.733 A in high-temperature 
flare spectra observed in RESIK channel 1, and are also available for analysis in the same 
way. The observations and analysis are described in Section 2, and comparison of argon 
abundances from the Ar XVII and Ar XVIII lines and with values obtained by other workers 
is described in Section 3. 



OBSERVATIONS 



The spectra discussed here were obtained in channels 1 (3.40-3.80 A) and 2 (3.83- 
4.27A) of the RESIK instrument. X-rays from the Sun were diffracted by silicon crystals 
(Si 111, 2d = 6.27 A) for both these channels. RESIK channel 2 spectra include Ar XVII 
emission lines consisting of three features called w (transition Is 2 1 S'o — ls2p 1 P\, 3.949 A), 
x + y (ls 21 S - ls2p 3 P 2)1 , 3.969 A), and z (ls 21 S - ls2s 3 S u 3.994 A). The line flux ratio 
G = (x + y + z)/w is slightly temperature-dependent. 

The G ratios from measured line fluxes in particular time intervals were plotted against 
single values of temperature and emission measure J v N 2 dV (N e = electron den sity V 



emitt ing volume) estimated from the ratio of emission in the two channels of GOES (IWhite et al. 

2005), and were found to be close to theoretical values based on data from chianti ( Sylwester et al. 



2008a|). Although a slight improvement to the fit was obtained by taking a two-temperature- 
component emission measure, in this work we take the original isothermal model as more 
suitable for analyzing the large number of spectra available. Channel 1 spectra include 
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K XVIII lines which were previously discussed (jSylwester et al.ll2010l ) and, for temperatures 
greater than about 10 MK, the Ar XVIII Ly-a line feature at 3.733 A. 

Some 2795 spectra from a sample of 20 flares between 2002 August and 2003 February, 
having GOES importance rangin g from CI to XI, were selected for anal ysis, and an average 
of non- flaring points taken from ISylwester. Sylwester. fc Phillips! (120101 ) . The spectra were 
collected over integration times that were inversely related to the X-ray activity level, from 
2 s at the peaks of s trong flares to 5 minutes for non - flaring per iods. Details o f the s elected 
spectra are given by ISylwester. Sylwester. fc Phillips! ( 120101 ) and lPhillips et al.l (j2010l ). Cali- 
bration factors, converting photon count rates to absolut e flux units, were applied to all the 
spectra from post-launch analysis (ISylwester et al.ll2005l ). 



The Ar XVII lines in channel 2 were observed in spectra over a large temperature range. 
As in our previous work and in flare analyses by several other authors, we derived temperature 
(Tgoes) from the flux ratio of t he two GOES chan nels for the time period of each RESIK 
spectrum using the procedure of I White et al.l (120051 ) . Figure [1] (left-hand plot) shows all the 
spectra analyzed stacked vertically in order of T GO es- in the right-hand plot, the spectra are 
normalized to the total number of photon counts, and the total spectrum shown at the top 
of each plot. The Ar XVII lines are evident over the entire temperature range shown, even 
at temperatures as low as 5 MK. As well as the principal Ar XVII lines (3.949-3.994 A), the 
S XV Is 2 1 Sq — lsAp 1 Pi line (also known as wA) is prominent at 4.088 A, as is a line feature 
at 4.19 A made up of S XIV dielectronic satellites (transitions Is 2 2/ — Is2l4l', I, I' = s, p) 
and at higher temperatures (T e > 10 MK) the H-like CI (CI XVII) Ly-a line. Weak lines of 
S XV with transitions Is 2 — lsnp with n ^ 10 occur between 3.998 A (wh) and 3.883 A (wlO) 
- wavelen gths are from atomic data extending the usual datab ase of the chianti database 
and code JDere et al.l[l~9~9~7l : JLandi et allboOfit bere elaH boOflh fl 

The Ar XVIII Ly-a line at 3.733 A falls in the range of RESIK channel 1. It is emitted 
over a wide range of temperatures, though it is most prominent a t temperatures above a bout 
10 MK. The line is visible in channel 1 spectra already shown in ISylwester et al.l ( 120101 ) in a 
stacked representation like Figure HJ The Ar XVIII line feature occurs at very nearly twice 
the wavelength of the well known group of highly ionized Fe lines (2 x 1.85 A) but these have 
no effect on the observed spectra as second-order diffraction is forbidden for Si 111 crystal 
used in channel 1. Weak S XVI lines (Is — 5p, 3.696 A; Is — 4p, 3.784 A) occur near the 
Ar XVIII line feature, but have no effect on the measurement of the Ar XVIII line flux. 



1 Thanks are due to Dr E. Landi for supplying these data. 
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3.9 4.0 4.1 % 4.2 
Wavelenght [A] 



3.9 4.0 4.1. 
Wavelenght [A] 



Fig. 1. — RESIK channel 2 spectra stacked with increasing temperature Tqoes ( m MK) 
in the vertical direction shown in a gray-scale representation (lighter areas indicate higher 
fluxes). The total spectrum is indicated at the top of each plot, with synthetic spectra 
(narrow line profiles) also shown. The principal lines are: Ar XVII w, x + y, z (3.95 A, 
3.97 A, 3.99 A); S XV wA (4.09 A), and the 4.19 A line feature, which is a blend of S XIV 
satellites and the CI XVII Ly-a lines. Left: Not normalized. Right: Normalized to the total 
number of photon counts in each spectrum. 
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ANALYSIS AND RESULTS 



The fluxes of the Ar XVII w and Ar XVIII Ly-a lines were estimated by taking narrow 
wavelength intervals around the lines (3.94-4.01 A for the w line and 3.72-3.745 A for 
the Ly-a line) and subtracting a portion of nearby continuum. Uncertainties in each flux 
measurement were taken from the uncertainties in the photon counts, assuming Poissonian 
statistics. For the Ar XVII w line, the continuum flux was invariably much less than the 
flux in the range containing the line, but this was not always the case for the Ar XVIII Ly-a 
line. We combined the uncertainties in the line and continuum flux measurements and if 
this combined uncertainty for any spectrum was larger than the difference between the line 
and continuum flux we rejected the estimate. 

To obtain estimates of the Ar abundance, we used a procedure similar to the one used 
by ISylwester et al.l ( 120101 ) for the abundance of K from K XVIII lines observed in channel 1 . 
GOES temperatures Tgoes an d emission measures EM GOES for the time interval correspond- 
ing to each spectrum were found and the value of the flux of each line divided by -EMgoes 
(units of 10 48 cm -3 ) was plotted against Tgoes- The Ar XVII w line plot is shown in Fig- 
ure |2] (left panel). Points from individual spectra are shown as black dots while averages over 
1-MK temperature intervals are also shown. The much smaller sca tter of points than for the 
corresponding plot for K XIX lines (Fig. 3 of ISylwester et al.l (120101 )) reflects the larger fluxes 
of the Ar XVII line. To these points are added theoretical contribution functions G(T e ) given 
by 



G(T e ) 



iV(Ar+ ib ) iV(Ar +lb ) iV(Ar) JV(H) 
iV(Ar+ 16 ) iV(Ar) JV(H) N e 



N P A 



;o 



cm s 



;i) 



where N(Arf 16 ) = population of the excited level i = ls2p 1 P 1 of ion Ar^ 1 " 16 , iV(Ar +16 ) = 
number density of the ion Ar +16 (all levels summed), iV(Ar) = the number density of all 
ionization stages of Ar, A(H) = number density of hydrogen (H), N e is the el ectron density, 
and A i0 is the transition probability from level i to the ground state. Following iPhillips et al. 



(120081 ). N(H)/N e was taken to be 0.83, which is correct for coronal plasmas with T e > 10 5 K. 
The atomic data, taken from versio n 6 of chianti, include excitation rate coefficients from 
the Ar +16 ground state, taken from IWhiteford et al.l (1200 ll ). based on -R-m atrix code calcu 



lations of Ar +16 up to n = 4 levels including effects of radiation damping. IWhiteford et al. 



(j200ll ) estimate the uncertainties in the rate coefficients for excit ation to the 1s2jp l Pi level 
to be ±10%. The ion fractions N (Ar +m ) / N (Ar) were taken from brvans et al.l (120091). We 



took the recent Ar a bundance estimates N (Ar)/N(H) = A(Ar) from both iLodderd ( 20081 ) 
(A(Ar) = 6.50) and Asplund et all (120091 ) (A(Ai) = 6.40): different line-styles show the 
G(T e ) curves for the two assumed abundances which are shown in Figure [21 The contri- 
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bution of unresolved Ar XVI dielectronic satellites, which is accounted for in chianti, is 
included in the G(T e ) curves. 

A similar procedure was used for flux estimates of the Ar XVIII Ly-a line, where again 
GOES temperatures and emission measures were used. The resulting plot is shown in Fig- 
ure [3] (left panel). The theoretical contribution functions were also derived from chianti 
and include unresolved satellite lines. In t his case, excitat i on rat e coe fficients are taken from 



i nterp olation of i?-matrix calculations of lAggarwal et al.l ( 119921 ) and lAggarwal fe Kingston 



( 119931 ) for ionized He and Fe. This is a rather wide stren gth of atomic number s, and there 



are now more recent data for specifically Ar XVIII lines (lAggarwal et al.l 120081 ) but as the 
transitions are simple, it is thought that interpolation for the case of Ar XVIII [Z = 18) will 
not lead to large errors. 

Each value of the observed line flux divided by EMqqes gives a value for the Ar abun- 
dance N(Ar)/N(il), so a number distribution can be plotted for the abundance values ob- 
tained. This is done in the right-hand panels of Figures [2] and The derived value of 
A(Ay) = 12 + log 10 [N (At) / N (H.)} is plotted along the horizontal axis, and the numbers of 
spectra with particular values of A(Ay) over intervals of 0.1 are plotted along the vertical 
axis. An approximately gaussian distribution is apparent for the Ar XVII w line results, 
with the peak of the best-fit gaussian at v4(Ar) = 6.44 and the FWHM points at 6.35 and 
6.52. The value of A(Ay) from the Ar XVII line fluxes is therefore 6.44 ± 0.07 (uncertainty 
is standard deviation). The Ar XVIII Ly-a line is formed at a higher temperature than the 
Ar XVII w line, and in fact the uncertainties in the wavelength interval defining the line 
(3.72-3.745 A) become larger than those in the nearby continuum for Tqoes < 15 MK. We 
thus took only line flux measurements for T GO es ^ 15 MK to derive the argon abundance. 
These estimates are shown in the right-hand panel of Figure [3j where as with the argon 
abundance estimates from the Ar XVII w line an approximately gaussian distribution is 
obtained. The peak value A(Ay) = 6.49, with FWHM points at 6.30 and 6.68. The value of 
A(At) from the Ar XVIII line fluxes is therefore 6.49 ± 0.16. The larger uncertainty reflects 
the weaker flux of the Ar XVIII Ly-a and the smaller number of flux measurements. The 
mean of the Ar XVII and Ar XVIII values weighted by the variances is A(Ar) = 6.45 ± 0.06. 

There is apparently little variation from flare to flare in the derived values of ^4(Ar) from 
the Ar XVII or Ar XVIII lines. This can be seen from Figure H] (top panel) where the mean 
value of A(At) is plotted for each of the 20 flares here, together w ith a non-flaring period 



discussed in previous work (jSylwester. Sylwester. fe Phillips! 120101 ) . The v alue of A(K) is 



also p lotted, being derived from the K XIX lines in channel 1, as reported by lSylwester et al. 



(J2010l ). as well as the logarithm of the abundance ratio, i.e. A(K) — A(Ar). The error bars 



indicate the uncertainties in the mean values. The total variations are within the estimated 



G(T) for Ar band 3.94-4.01 A 



A Ar histogram 




15 
T GO es [MK] 



6.0 6.2 6.4 6.6 6.8 7.0 
Log Abund [A H =12] 



Fig. 2. — (Left panel:) Estimated fluxes (small black points) in the Ar XVII w line divided 
by the GOES emission measure EMqoes (units of 10 48 cm~ 3 ) plotted against Tgoes (an 
isothermal plasma is assumed). Small circles (yellow) are means of the points over 1 MK 
intervals. The dotted and dashed curves are the theoreti cal G(T P ) curve s for the Ar XVII w 



l ine ba sed on chianti with the Ar abundanc e estimates of iLodderd (120081 ) and lAsplund et al. 



(120091 ) respectively, and the ion fractions of iBryans et al.l (120091 ) . They include unresolved 
satellite lines. (Right panel:) Histogram of arg on a bundances (A(Ar)) in intervals of A(Ar) = 
0.1, with the Ar abundances of lLodderd (120081 ) and lAsplund et al.l (120091 ) indicated as vertical 
lines (dotted and dashed respectively). The peak of the distribution gives A(Ai) = 6.44, 
and the distribution width (full width half maximum) d efines th e unce rtainty range (6.35- 
6.52). Blue curves an d vertical lines refer to the work of iLodderd (120081 ). red to the work of 
Asplund et all (l2009h . 
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Fig. 3. — (Left panel:) Estimated fluxes (points) in the Ar XVIII Ly-a line feature at 3.74 A 
divided by the GOES emission measure EM GOES (units of 10 48 cm -3 ) plotted against T G qes- 
As with the analysis of the Ar XVII lines, an isothermal plasma is assumed for al l spectra. 
The theoretical curyes are from chianti with the Ar abundances of lLodderd ( 20081 ) (dotted) 
and lAsplund et al.l ( 120091 ) (dashed). (Right panel:) Histogram of argon abundances (A(Ar)) 
plotted as in Figure for spectra having T G qes > 15 MK. Here the distribution gives 
A(Ar) = 6.49 with an uncertainty range from the distr ibution width (FWHM) of 6.30 to 
6.68. Vertical l ines i ndicate the argon abundance from lLoddersI ( 120081 ) (dotted) and from 
Asplund et al. (120091) (dashed). Blu e curves and vertical lines refer to the work of lLodders 
fl2008r ). red to LAsplund et all (J2009h . 
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uncertainties, and so the abundances of K and Ar appear to be constant fo r this sample of 
RESIK flares. The non-flaring points from ISylwester. Sylwester. &: Phillips! (120101 ) are con- 
sistent with the remaining flare points though the uncertainties are somewhat larger because 
of the weaker K XIX and Ar XVII line emission at lower, non-flare temperatures. Possible 
time variations of the K and Ar abundances were checked because previous results from So- 
lar Maximum Mission Bent Crystal Spectr ometer data ind i cated time- varying abundances 
of Ca from the Ca XIX w line (3.177 A: see ISvlwester et~aD (jl984h )- 



4. DISCUSSION AND CONCLUSIONS 



Our values of the Ar abundance A(Ar) from RESIK Ar XVII w line, 6.44±0.07, and from 
the Ar XVIII Ly-a line, 6.49 ± 0.16, are very close to each other, the small difference being 
much less than the uncertainties indicated by the scatter of individual points in Figures |2] and 
EJ The uncertainties in the excitation rate coefficients for the Ar XVII lines are, as mentioned, 
±10%, which translates to a vertical shift of the theoretical contribution function G(T e ) by 
this amount. This is very small compared with the scatter of measured values, however, 
even for the more intense Ar XVII w line. Even so, a slightly better estimate m ight result 
from use of the Ar XVII i?-matrix rate coeffi cients of lAggarwal fc Keenanl (J2005[ ) and more 
particularly the Ar XVIII rate coefficients of lAggarwal et al.l ( 120081 ) since they are likely to 
be an improvement over the interpolated values presently used in the chianti database for 
this ion. 

Our abundance est imates may be comp ared with other published values. The estimate 
from RESIK spectra of iPhillips et al.l (120031 ) is based on spectra from the 2002 period when 
the instrument characteristics were not quite optimized. In addition, a co nstant intensity 
calibr ation factor for each channel was taken though the more recent work of ISylwester et al. 
(200 5J) shows that there is a slight variation with wavelength. However, the value from 
Phillips et aD (l2003h . A(Ax) = 6.45±0.03, is very close to the value obtained here (6.44±0.07) 
from the Ar XVII w line. 

Our values are somewhat less (up to 0.20 dex) than values ob tained from Skylab ul- 
traviolet lines ( JFeldman fc Widingl Il992l ; iFeldman fc Laming! |2000|) but are very close to 
that obtained from P78-1 observations of the Ar XVII Is 2 — ls5p X-ray line at 3.128 A 



( Doschek et al.lll985l ). The 3.128 A line used was not generally visible with the PI 8-1 spec- 
trometer but became apparent when t he instrument was sl ightly off-pointed; the quality 
of the spectra are reasonably good and iDoschek et al.l ( 119851 ) consider their estimate to be 
reliable. The SMM Flat Crystal Spectrometer scanned through the wavelength region of the 
Ar XVII lines at 3.93 A during flares on a few occasions; the spectra, which do not have 
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Fig. 4. — Estimated abundances of K, A(K), from K XVIII line w (3.53 A) and Ar, A(Ar), 
from Ar XVII line w (3.95 A), and the logarithm of the K/Ar abundance ratio, A(K)—A(Ar). 
For comparison an average of non- flare points is added (last point in each plot). 
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particularly good statistical quality, are discussed by [Phillips et al.l (J199J). Observed ratios 
of the Ar XVII w line with the nearby S XV Is 2 — ls3p (w3) and Is 2 — ls4p (w4) lines gave a 
value of the Ar/S abundance ratio of 0.3 ± 0.1, which with the coronal value of A(S) = 6.99 
from ISylwester et al.l (J2008bl ) leads to A(Ar) = 6.47 ± 0.15, in agreement with the values 
obtained in this work. 



The recent compilations of abundance data of lLodderd (120081 ) and lAsplund et al.l (120091 ) 
give A(At) = 6.50 ± 0.1 and A(Ar) = 6.40 ± 0.13 respectively, the chief difference being the 
abundance of O taken in these two works. Our value from Ar XVII w line flux measurements 
falls between these two values and agrees with both to within the stated uncertainties. 

Some discussion on the nature of the FIP effect has resulted in different conclusions 



about coronal and proxies for photospheric argon abundances: iFeldman et al.l (119921 ) state 
that high-FIP (FI P > 10 eV) elements have approximately the same abundances whereas 
Fludra fc Schmela (119991 ) argue that the coronal abundance of high-FIP elements are de- 
pleted by a factor 2. Considerably more data on ar gon abun d ances are now available since 
either of these publications, and are summarized by lLoddersI (120081 ) . A mean of her values 
of the proto-solar abundances, i.e. those assumed to be appropriate to the solar nebula, 
from a variety of objects (H II regions, Jupiter's atmosphere, planetary nebul ae, and B 
stars) is A(Ax) = 6.55 with a scatter of about ±0.1. Following lLoddersI (120081 ). reducing 
this abundance by 0.07 to allow for heavy-element settling over the Sun's lifetime gives a 
solar photospheric abundance of 6.48 which again is in agreement with our values. This 
suggests that there is little or no depletion of the coronal Ar abundance with respect to 
the photospheric, assuming that the adjusted proto-solar abundances are a reflection of the 
latter. 
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